Parity anomaly and Landau-level lasing 
in strained photonic honeycomb lattices 
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We describe the formation of highly degenerate, Landau-level-like amplified states in a strained 
photonic honeycomb lattice in which amplification breaks the sublattice symmetry. As a consequence 
of the parity anomaly, the zeroth Landau level is localized on a single sublattice and possesses an 
enhanced or reduced amplification rate. The selection of the sublattice depends on the strain orien- 
tation, but is independent of the valley; therefore, the anomaly directly affects the mode selection 
and lasing threshold. The spectral properties of the higher Landau levels are constrained by a 
generalized time-reversal symmetry of dynamical origin. 

PACS numbers: 42.55.Tv, 03.65. Vf, ll.30.Er, 73.22.Pr 



Nonuniform deformations of the honeycomb lattice of 
graphene result in a pseudomagnetic field which deflects 
particles in analogy to the Lorentz force, with small 
amounts of strain producing fields that are large enough 
to create well-defined Landau levels in the low-energy 
range of the spectrum [D-Q- Here we describe how the 
addition of gain in an analogous photonic setting results 
in the formation of highly degenerate amplifying Landau 
levels, which can provide the platform for a laser with 
macroscopic mode competition. 

The spectral properties of these levels become intrigu- 
ing when the gain breaks the sublattice symmetry. Be- 
cause of the parity anomaly 0^3; the amplification of 
the zeroth Landau level is dictated by one of the two sub- 
lattices, which here is selected depending on the strain 
orientation. A reflection symmetry entails that the in- 
stances of this level in the two k-space valleys are de- 
generate. In contrast, the time-reversal symmetry T and 
the inversion symmetry V are broken by the amplifica- 
tion imbalance; these symmetries usually serve to can- 
cel the anomaly when both valleys are taken into ac- 
count ff| 0] ■ We also identify a generalized time-reversal 
symmetry which constrains the spectral properties of the 
higher Landau levels. The amplification rate of these lev- 
els equals the average rate of the two sublattices, up to a 
finite threshold of the imbalance at which two levels coa- 
lesce and their rates bifurcate. The consequences for the 
lasing threshold of the system are outlined at the end of 
this work. In particular, we find that depending on the 
strain orientation and signature of the amplification im- 
balance, the lasing threshold is set either by the zeroth 
or by the first Landau level. 

Model of a strained active photonic honeycomb 
lattice. — We specifically consider the active photonic sys- 
tem sketched in Fig. Q] Panel (a) shows a segment of a 
honeycomb lattice, with vertices which could represent 
weakly coupled optical fibers in a three-dimensional pho- 
tonic lattice [H, or, more abstractly, a set of basis states 
for a suitable band in a two-dimensional photonic crystal 
0, flol | ■ The honeycomb lattice consists of two sublat- 
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FIG. 1. (Color online) (a) Segment of a honeycomb lat- 
tice, with vertices representing states in a two-dimensional 
photonic crystal or weakly coupled optical fibers in a three- 
dimensional setting, (b) Sketch of a deformed arrangement 
which results in a constant pseudomagnetic field, (c) We in- 
vestigate the interplay of this field with amplification and 
absorption that breaks the sublattice symmetry. The two 
sublattices A and B have amplification rates 7 a and 7b, re- 
spectively (negative values correspond to absorption). The 
pseudomagnetic field resulting from the strain is modeled via 
smooth coupling functions ti whose definition ([5| involves the 
bond vectors p l , I = 1, 2, 3. 



tices, A sites and B sites, which we equip with different 
amplification or absorption rates, in analogy to recent 

- lH of non-hermitian 
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works on optical realizations 

■PT-symmetric quantum mechanics~jl5| . (In the present 
setting, V stands for the inversion about the center of a 
hexagon, which maps A sites to B sites and thus inverts 
the amplification imbalance; T corresponds to complex 
conjugation and converts amplification into absorption, 
which also inverts the imbalance.) Panel (b) sketches 
an inversion-symmetry-breaking deformed arrangement 
which results in a constant pseudomagnetic field whose 
interplay with the symmetry-breaking effects of ampli- 
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FIG. 2. (Color online) Dependence of the Landau level spec- 
trum ((9]), (|10[) on the amplification imbalance 7 = (74 — 
7s)/2, for strain leading to a pseudomagnetic field of strength 
p. (a) Real part, which vanishes for the zeroth Landau level, 
as well as for the other Landau levels beyond their bifurcation 
thresholds 7„, Eq. (b) Imaginary part, which becomes 

finite beyond the bifurcation. Because of the parity anomaly, 
the zeroth Landau level breaks the symmetry of the spectrum 
as it is located on the A sublattice for /3 > 0, resulting in the 
thick solid line, while it is located on the B sublattice for ft < 
(thin dashed line). The circles indicate the lasing threshold 
when I7I is increased at fixed average absorption 7 = —71. 
For 7 > this threshold is set by the zeroth Landau level, 
7l = 71 (solid circle), while for 7 < it is set by the first 
Landau level, jl = —V^Ji (open circle). 



((Pa(t), 1 Pb(^)) T which is obtained by separating out 
rapid fluctuations with wave vector K CT = <T(47r/3v3p)i, 
a = ±1, where a distinguishes two independent valleys 
whose existence is guaranteed by the V and T symmetries 
of the unstrained passive system [10(. The eigenvalues e 
of H determine the frequencies u) of quasibound states in 
a two-dimensional setting [3, [13] or the propagation con- 
stant ck z along the third direction in a three-dimensional 
setting Eigenvalues with a positive imaginary part 
correspond to amplified states (in time or along the prop- 
agation direction) , while those with a negative imaginary 
part correspond to decaying states. 

For vanishing 7a = 7s = and constant t a b = to, the 
system is periodic and the band structure displays the 
familiar Dirac cones e = ±u|q| near each corner of the 
Brillouin zone (the K and K' points situated at K + and 
K_, respectively), where q = k — K ff is the wave vector 
relative to the corner point [16| . Weak uniform strain, 
with t a b — £1,2,3 only depending on the bond orientation, 
displaces the cones from the corners by an amount A 
[E 1|, E3 • In the presence of amplification and absorp- 
tion with = — 1b = 7, the full band structure of the 
uniformly strained system can still be real since in this 
case the non-hermitian Hamiltonian ([1} displays the VT 
symmetry 13[ 14 j 



fication and absorption we are interested in. Panel (c) 
illustrates the microscopic modeling of these effects. The 
sublattices carry amplification rates ^a — 7 + 7 an d 
7b =7 — 7, respectively, where 7 is the average rate 
and 7 quantifies the imbalance. The rates ^a and jb 
may be negative, in which case they signify absorption. 
Strain results in a spatial variation of the coupling terms 
tab between neighboring A and B sites a and 6, respec- 
tively, which we parameterize as t a b — ti(r a ), where r a 
is the unstrained position of the A site and I = 1,2,3 
indicates the orientation along the unstrained bond vec- 
tors Pi (with \pi\ = p the unstrained nearest- neighbor 
distance). The typical magnitude of the coupling terms 
is denoted as to- 

We focus on a spectral range where the unstrained pas- 
sive lattice displays a conical band structure, mimicking 
graphene [13, EB] ■ Jhi combination of the descriptions of 



strained graphene 
tices or crystals [1. 



li. 211171 a nd photonic honeycomb lat- 
10 . ll3L Il4j , the resulting Lorentz force 
and the effects of amplification and absorption are then 
captured by a Dirac equation with Hamiltonian 



H = 



11A 

v{aP x +iP y ) 



v(aP x - iP y 
«7b 



(1) 



v = Mop/2, P x = -id x - A x , P y = -id y - A y , where 
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(*2 - t 3 )j (2) 



is the pseudomagnetic vector potential. This Hamil- 
tonian applies to a continuous spinor wave function 



H(x,y) = a x H* 



-y)<j x = VTH{x,y)VT, (3) 



where a x is the Pauli matrix. The VT symmetry is 
spontaneously broken when 7 exceeds a threshold where 
eigenvalues coalesce and the associated eigenstates cease 
to be joint eigenstates of VT; this leads to complex 
branches of the band structure where e(k) is complex 

3 El- if 

amplification and absorption are imbalanced, 
all eigenvalues are shifted by V). This includes the case 
of 'passive' VT symmetry, where 7 = — I7I such that one 
sublattice is absorbing and the other sublattice is neu- 
tral. In these more general slightly relaxed VT 
symmetry can be stated as 



H = VTUVT + 2«7. 



(4) 



The spectrum of such a system is constrained to eigen- 
values which either fulfill Ime„ = 7, or are paired with 
another eigenvalue = e n — 2ilme„ + 2«7. However, 
strain explicitly breaks the VT symmetry, as we explore 
in the following. 

Landau levels. — We consider a strain configuration 
which results in a constant pseudomagnetic field of 
strength f3. This follows from a smoothly varying three- 
fold symmetric configuration with [l[ 

t,=to[l-G8/2)p,.r], 1 = 1,2,3, (5) 

which gives rise to a vector potential 

A = (aP/2)(-yi + xj). (6) 
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Microscopically j3 depends on the sensitivity of the cou- 
pling terms on the nearest-neighbor spacing, as well as 
on the strain orientation; here we assume that this pa- 
rameter is given. Assuming unless otherwise stated that 
P > we can write the Hamiltonian as 

v^zpn i lB J' [) 

1 - r (-iapx/2 - iad x + py/2 + d y ), (8) 



u 



n : 



2/r 



where [II, 11^] = 1 coincides with the algebra of harmonic 
oscillator annihilation and creation operators (the same 
algebra applies to both valleys, a = ±1, which becomes 
important below). This delivers a spectrum of Landau 
levels, with the zeroth level given by 



Xo 


a/2 



£q = ilA = «7 + *7) 



(9) 



where xo = {P/2tt) 1/2 exp[-/3(ar + y )/4 + \(ax + iy) - 
A /P\ represents the infinitely degenerate set of Landau 
states fulfilling ILxo = 0. With X m = (m!)" 1 ' 2 ^)^, 
m = 1, 2, 3, . . ., the other Landau levels are given by 



X\n\ 
Oi n X\n\-l 



i7 + sgn(ri) v / 7 I 2 l -7 2 , (10) 



sgn(nWl-^-^ 5 7 n = V2v 2 P\n\, (11) 
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for n = ±1, ±2, ±3, . . .. This spectrum is shown in Fig. [2] 
At vanishing ja,b (thus 7 = 7 = 0), these solutions 
reduce to the strain-induced Landau levels that have 
been studied in the graphene literature [1-4]. For uni- 
form amplification or absorption (7 = 0), the levels are 
simply shifted by 17. For a finite amplification imbal- 
ance (7^0), the levels with opposite level index n and 
n = —n (with n 7^ 0) coalesce at a threshold value 
I7I = 7„ and then bifurcate into a pair of levels which 
fulfill the spectral constraints stipulated below Eq. (QJ; 
in particular, the average imaginary part of these eigen- 
values is given by i"f. The zeroth Landau level, however, 
has an imaginary part which differs from 7, and is not 
accompanied by a partner state (not even in the other 
valley). This feature of the zeroth Landau level rules out 
the existence of any 'PT-like antiunitary operator which 
would commute with the Hamiltonian. 

The special nature of the zeroth Landau level can be 
seen as a direct consequence of the parity anomaly [H-0] , 
which in the present context is most conveniently identi- 
fied by considering the supersymmetric interpretation of 
Hamiltonians of the form (JT]) 5] . Depending on whether 
one eliminates the B site or A site wave function, re- 
spectively, the corresponding eigenvalue equation can be 
written as 

(e - i-ia)va = (e - nsy^Pv^n^A, (12a) 
{e-i lB ) VB = (e-i7A)~ 1 2/?w 2 nnVs, (12b) 



which provides a simple example of supersymmetric part- 
ner potentials (since ITlL = IFII + 1, the potentials 
are simply shifted with respect to each other). Both 
equations deliver the same spectrum, except for the ze- 
roth Landau level, which only occurs in the spectrum 
of Eq. (|12ap as ipB = is not normalizable. This state 
without an iso-spectral partner thus breaks the sublattice 
symmetry — its wavefunction is localized on the A sub- 
lattice, and for 7^4 ^ jb this asymmetry is reflected by 
a departure from the overall symmetry of the spectrum 
about £7. This holds for P > 0, as we have assumed 
so far. For P < 0, one needs to modify the definitions 
of II and nt such that in the Hamiltonian ([7]) they are 
effectively interchanged, and the zeroth Landau level is 
localized on the B sublattice, with eo = i"fB- 

Focussing on a single valley (say around the K point, 
a = 1), this anomaly is fully analogous to the parity 
anomaly in the problem of massive Dirac electrons in a 
magnetic field, which possess an extra state located at en- 
ergy E — mc 2 or E = —mc 2 (depending on the sign of the 
field) that breaks the symmetry of the spectrum about 
E = [7]. For electrons on an ordinary honeycomb lat- 
tice this anomaly is canceled in the K' point [a, 0] , which 
can be related to the K point by either using the V sym- 
metry (which interchanges the two sublattices) or the T 
symmetry (which inverts the magnetic field). An uncom- 
pensated anomaly can only occur when these symmetries 
are broken, as is the case in the Haldane model, which 
includes suitable magnetic fluxes 0] • In the present pho- 
tonic setting the T operation relating the two valleys in 
k-space inverts the sign of the amplification imbalance 
7; furthermore, the V operation not only inverts 7 but 
also the direction of the vector potential Q — thus, both 
symmetries are indeed broken. Instead, the Hamiltonian 
([7| can be mapped from one valley to the other by the 
reflection symmetry x — > —x, a — > —a. Therefore, the 
parity anomaly for the zeroth Landau level is replicated 
identically in both valleys. 

We now turn to the other Landau levels. These 
are constrained by the particle-hole like symmetry 
H(x,—y)* = —1-L{x,y), which results in the pairing 
e_„ = — £* of eigenvalues before the bifurcation thresh- 
old, I7I < 7„ (after the bifurcation, | — / 1 > 7„, no con- 
straint arises from this symmetry since e n = — e* n when 
eigenvalues are imaginary). The question now arises why 
these levels also obey the spectral constraints that are 
usually associated to an antiunitary T'T-like operator 
commuting with a non-hermitian Hamiltonian. These 
constraints imply Ime„ = Ime_„ = 7 before the bifur- 
cation, while after the bifurcation the eigenvalues come 
to lie symmetrically with respect to the line Ime = 7. 
Furthermore, we observe that before the bifurcation, the 
associated wave functions (fT0|) each have equal weight on 
the A and B sublattices, |a n | = |a_„| = 1, even though 
the spatial dependence oc X|n|iX|n|-i is different. Af- 
ter the bifurcation, |a„| = l/|a_„| ^ 1, and the level 
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with the larger imaginary part has a larger weight on the 
more amplifying sublattice, while the other state is pre- 
dominantly localized on the opposite sublattice. These 
properties are all compatible with^ existence of a gen- 
eralized time-reversal symmetry VT applying to these 
states. Such a symmetry needs to meet the following de- 
mands: VT must commute with the Hamiltonian in the 
sense of Eq. (jH), it must square to the identity, and it 
must act as an antiunitary operator in the subspace of 
the higher Landau levels, while not mixing these with the 
zeroth Landau level. 

To identify this symmetry we introduce the basis 

\m,A)=( x ™\ \ m> B)=(^y m = 0,1,2,... 

of ordinary Landau states localized on the A or B sub- 
lattice (we suppress the degeneracy of these levels). In 
this basis, the Hamiltonian ([7]) takes the form H = 
ij A \0,A)(0, A\ +H, where 

OO 

U = (hA\m + l,A){m + l,A\ + i"/ B \m,B}{m,B\ 

m=0 

+ vy/2P{m + l)(|m + 1, A)(m, B\ + \m, B){m + 1, A\)j 

is the Hamiltonian in the subspace excluding the zeroth 
Landau level |0, A). Inspecting the properties of 

OO 

V = ^T(\m + l,A){m,B\ + \m,B)(m + 1, A\) (13) 

m=0 

f:T\m,L)^T*\m,L), L = A,B, (14) 

in the original Hilbert space, we find T 2 = 1, P> = V, 
V 2 = 1 - |0,A)(0,A|. Thus VT = Vf\ n ^o is an an- 
tiunitary operator in the space of higher Landau lev- 
els. Furthermore, an explicit cakvulation now delivers 
the desired relation H = VTHVT + 2^7 in the space 
of these levels, which entails the spectral constraints. 
This symmetry is of dynamical origin as its construc- 
tion makes explicit reference to the eigenstates of the 
system. The symmetry does not extend to the zeroth 
Landau level since V is not unitary if this state is in- 
cluded. _ While VTH = HVf - 2i^vf, the relation 
VTUVT = H-2ij- 2ij\0,A)(0,A\ again reveals the 
special spectral status of this level. 

Conclusions. — In this work, we mainly focussed on the 
symmetry aspects that govern the spectrum of Landau 
levels in the described strained photonic lattice with am- 
plification that breaks the sublattice symmetry. One can 
now consider the onset of lasing, which occurs when one 
of the Landau levels acquires a positive imaginary part. 

For fixed 7 < 0, the system is passive at 7 = (uniform 
absorption), but as I7I increases the zeroth Landau level 
changes its imaginary part, and so do the other Landau 
levels beyond their bifurcation thresholds j n . The lasing 



threshold 7^ now depends on the sign of 7. If 7 > 0, 
the lasing threshold is given by 7^ = I7I since the zeroth 
Landau level then is located on the amplifying sublat- 
tice. For 7 < 0, the first level to meet the real axis is 
associated with the pair n = ±1 involving the first Lan- 
dau level, with lasing threshold 7l = \/2v 2 (3 + 7 2 [see 
Eq. (jlQI) ]. In the case 7 = 0, the lasing threshold either 
vanishes (for 7 > 0) or is finite (for 7 < 0), depend- 
ing on whether it is set by the zeroth or first Landau 
level. These considerations apply to the principal strain 
orientation studied here ((3 > 0). If the (3 is set to a nega- 
tive value, the role of the two sublattices is interchanged 
and the zeroth Landau level becomes lasing for 7 < 0. 
Similar asymmetric threshold scenarios arise when one 
approaches lasing by changing 7 at fixed 7, or when one 
changes the amplification rate on one sublattice only. 

In the lasing regime, a macroscopic number of modes 
in the zeroth or first Landau level will participate in the 
mode competition. In a finite system, the exact degen- 
eracy will be lifted by the boundary conditions, but this 
lifting will be small in the bulk, while edge states can 
also appear. We leave the investigation of these aspects 
for future work. 

Note added. — Since the submission of this article, pseu- 
domagnetic fields in strained photonic systems have been 
realized experimentally in Ref. [lH; this work observes 
signatures of Landau levels and also suggests the possi- 
bility of degenerate lasing. 
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